Self-trapping transition for nonlinear impurities embedded in a Cayley tree 
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The self-trapping transition due to a single and a dirtier nonlinear impurity embedded in a 
Cayley tree is studied. In particular, the effect of a perfectly nonlinear Cayley tree is considered. 
A sharp self-trapping transition is observed in each case. It is also observed that the transition 
is much sharper compared to the case of one-dimensional lattices. For each system, the critical 
values of \ f° r the self-trapping transitions are found to obey a power-law behavior as a function 
of the connectivity K of the Cayley tree. 
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I. INTRODUCTION 

The interaction of an electron or an exciton with 
the lattice vibrations is of fundamental importance in 
understanding the electric properties of solids. For ex- 
ample, the transport of quasi-particles such as elec- 
trons or excitons in solids is highly influenced by the 
electron-phonon interactions. The consequences have 
been investigated using different methods Recently, 
these systems have been studied based on the rigor- 
ous analytical treatments and numerical solutions of 
simple models such as nonlinear Schrodinger equations 
|^|,D . One of them with varieties of applications in dif- 
ferent areas of science is the one-dimensional discrete 
nonlinear Schrodinger equation, given as || [l8| 



.dC n 
I-—— 
dt 



= V{C n+x + C n _i) + (e„ - Xn\C n \ 2 )C n , (1) 



where e n is the static site energy at site n and \n is 
the nonlinearity parameter associated with the nth grid 
point. Since J2 n I C n | 2 is set to be unity by choosing 
appropriate initial conditions, | C n | 2 can be consid- 
ered as a probability of finding a particle at the nth 
grid point. One way to derive this set of equations 
is to couple in the adiabatic approximation (in which 
the lattice oscillations are much faster than the exciton 
motion) the vibration of masses at lattice points of a 
lattice of N sites to the motion of a quasi-particle in 
the same lattice. The motion of a quasi-particle is de- 
scribed, however, in the frame work of a tight binding 
Hamiltonian (TBH). fn other physical context, the set 
of equations are often called the discrete self-trapping 
equations (DST). 

The analytical solutions of Eq. (1) are, in general, 
not known. However, for nonlinear quantum dimers 
which are two-site systems with the nonlinearity ei- 
ther on both the site-energies or in one of them can be 
solved analytically for any arbitrary initial condition. 
^From the analytical solutions, a self-tapping transi- 
tion is found in this model |^-[lC[| . The self-trapping 
transition for the symmetric dimcr is found at x/^=4 
||. The trapping of hydrogen ions surrounding oxy- 
gen atoms in metal hydrides and the energy transport 
from the absorption center to the reaction center in 
photosynthetic unit have been modeled by the effective 
quantum nonlinear dimer [|fj|- ^|JL^jr^ , [r7|] . The non- 
linear dimer analysis has also been applied to several 



experimental situations like neutron scattering of hy- 
drogen atoms trapped at the impurity sites in metals 
0, fluorescence depolarization M, muon spin relax- 
ation |^o| , nonlinear optical response of superlattices 
[ ^fj[ , etc. The self-trapping transition also occurs in 
the extended nonlinear systems rafiffl. An interesting 
experimental example in this context is the observation 
that trapped hydrogen atoms in metals like Nb move 
among the sites in the neighborhood of impurity atoms 
such as oxygen |20[] . All these studies have been per- 
formed for finite number of nonlinear sites by assuming 
that the quasi-particle is localized within the nonlinear 
sites. 

Later, the effect of nonlinear sites embedded in a 
host lattice on the dynamics of quasi-particles has been 
studied because of its importance in real systems. Dun- 
lap et al. f2l"|| studied the self-trapping transition due 
to a single nonlinear impurity embedded in an infinite 
one-, two- and three-dimensional host lattices. Self- 
trapping transitions were found at x/V=3.2, 6.72 and 
9.24 for one-, two- and three-dimensional simple host 
lattices respectively. Furthermore, the effect of the 
presence of a nonlinear cluster on the self-trapping 
transition has also been considered in one-dimensional 
host lattices Q. The study has also been extended 
to the case where the inertial effect of the lattice oscil- 
lators has been taken into account and rich trapping- 
detrapping transitions depending on the masses of os- 
cillators have been observed . All these studies have 
been made in one-dimensional host lattices. However, 
one needs to know whether or not the self-trapping 
transition occurs in hosts of different geometrical struc- 
ture. It would also be interesting to note the differences 
compared to the results for one-dimensional systems. 

The Cayley tree is one possible example of host lat- 
tices with different geometrical structure. An impor- 
tant variable characterizing the geometry of the Cay- 
ley tree is the connectivity K, which is the number one 
smaller than the coordination number, i.e., K = Z — 1, 
Z being the coordination number. The Cayley tree re- 
duces to a one-dimensional chain if K = 1. The struc- 
ture of the Cayley tree will be described in the next 
section. 

In this work, we study the self-trapping of an exci- 
ton in the Cayley tree with a single as well as a dimer 
impurity embedded in it. We also consider the fully 
nonlinear Cayley tree to observe the self-trapping ef- 
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feet. 



II. FORMALISM 

The structure of a Cayley tree with connectivity K = 
2 is shown in Fig. 1. For K = 1, the system reduces to 
a one-dimensional chain. In case of a single nonlinear 
impurity embedded in the Cayley tree, the neighboring 
sites surrounding the impurity site are symmetric. We 
thus consider the symmetric shells around the impurity 
site. The impurity site is designated as the zeroth site, 
and the nth shell is denoted by n, where n takes the 
values of 1, 2, 3, • • •, as the sites are away from the 
zeroth site (see Fig. 1). 




FIG. 1. The Cayley tree of the connectivity K = 2. The 
impurity is embedded at the site marked by 0. The sites 
in the first shell are numbered as 1, the sites in the second 
shell as 2, and so on. 

Thus, the nth shell contains ZK n ~ Y sites. All sites 
in the lattice have three nearest neighbors. While all 
the nearest-neighbor sites of the zeroth site fall in the 
first shell, two of the nearest-neighbor sites of any site 
in the nth shell fall in the (n + l)th shell and the rest 
one falls in the (n — l)th shell. We further notice that 
all sites in the same shell have an equal probability am- 
plitude. Therefore, under the tight-binding formalism, 
the time evolution of a particle (initially placed at the 
impurity site) on the Cayley tree may be governed by 
the following equations: 

i^L = ZC 1 - X \C \ 2 C 
dC 

= KC n+1 + n>l (2) 

where Co is the probability amplitude at the zeroth 
site and C n for n > 1 represents the probability am- 
plitude at any site in the nth shell. The x represents 



the nonlinear strength at the zeroth site of the Cayley 
tree. Without loss of generality, we take the nearest- 
neighbor hopping element to be unity. The normaliza- 
tion condition for the site amplitudes is given by 

oo 

\C \ 2 + Y,ZK"~ 1 \Cn\ 2 = l. (3) 

71=1 

Therefore, to observe the self-trapping transition due 
to a single nonlinear impurity, Eq. (1) should be solved. 

For the perfectly nonlinear Cayley tree, the time evo- 
lution for the site amplitudes on the Cayley tree of a 
particle (initially placed at the zeroth site) may be gov- 
erned by 

dt 
dC 

i__ = KCn+1 + Cn^ ~ x\C n \ 2 C n: n > 1 (4) 
at 

because the system may be treated in a similar way as 
it is treated for the single impurity case. 




FIG. 2. The Cayley tree with the connectivity K = 2. 
The dimeric impurity is embedded at the sites marked by 
and 1. 

If there is a dimeric impurity in the system, the sym- 
metry about one of the impurity sites does not hold 
anymore. The system, however, remains symmetric 
about the bond connecting the two impurity sites. In 
this case, the Cayley tree with dimeric impurity may be 
transformed to a one-dimensional system which might 
be studied more conveniently. The transformation has 
been reported in the earlier work [ p4[ ; however, for 
completeness, we briefly describe it in what follows. 

We pick a bond and assign the numbers and 1 on 
its two ends. The neighboring sites of the site 1 are 
numbered with an increasing order and those of the 
site are numbered with a decreasing order, as shown 
in Fig. 2. Thus, all points with the same number are 
assumed to be in the same generation and, accordingly, 
the number of points in the nth generation is K"^ 1 if 
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n > 1 and K> n > if n < 0. We note that all sites in a 
given generation have the same probability amplitude. 

We now consider the motion of a particle on a Cay- 
ley tree of the connectivity K with a dimeric impurity 
embedded at sites and 1. In the tight- binding for- 
malism with nearest-neighbor hopping only, equations 
governing the motion of a particle are 



.dC n 
I—— 

dt 

.dC n 
I—— 

dt 
.dCi 

dC 
i—— 

dt 



KCn+l+Cn-l, n>l 

A"C_|„|_i + C_|„|+i, n < 
KC 2 + C + e 1 C 1 
KC-i +Ci+ e Co, 



(5) 



where C n denotes the probability amplitude at any 
point in the nth generation. We assume that the ma- 
trix element of the nearest-neighbor hopping is unity, 
and that all points in a given generation arising due to 
a specific organization have the same site energy. The 
normalization condition for the site amplitudes gives 



oo 



(6) 



71=1 



We now make the following transformations: (i) r = 
VKt, (ii) eo = e /\/K and e x = ei/y/K, (iii) C n = 
K -(n-i)/2 Cn f 0r n > i an d C_ H = K~\ n \/ 2 C n for 
n < 0. Substituting these in Eq. (5) we obtain 



.dC„ 



i— — = C n+ \ + C„_i, for n > 1 and n < 0. 
dr 

■ dC\ ~ 1 ~ ~ 
i— = C 2 + -f=C + ad, 

dr JK 



■ dCa ~ 1 ~ „ ~ 

i-r- = C_i + -=Cx + e C . 
dr y/K 



(7) 



^From Eq. (7), the normalization condition reduces 

to X)^°oo |C"| 2 = 1- Therefore, the motion of a par- 
ticle on a Cayley tree is mapped to that on a one- 
dimensional chain. However, in this chain, the nearest- 
neighbor hopping matrix element between the zeroth 
and first site is reduced from unity to \j\f~K. It can 
be shown that the Green's function Go,o(E) calculated 
from Eq. (7) would yield Gq,o(E = E\J~K) for a Cayley 
tree of the connectivity K. Here, in the dimeric impu- 
rity problem, the impurities are defined as To = x|Co| 2 
and Tl = x|Ci| 2 with \ = ;t^- 

Eqs. (2), (4) and (7) can not be solved analyti- 
cally and, therefore, the numerical method, namely, 
the fourth order Runge Kutta method is employed. 
Since there is a conserved quantity in each case, the 
normalization condition is checked at every step of our 
numerical calculation. The time interval St = 0.001 
is used during the calculation. There are two ways to 
observe the self-trapping transition; One way is to look 
at the behavior of |C„| 2 as a function of t for various 



values of the nonlinear strength, and another is to look 
at the time averaged probability of the particle at site 
n as a function of the nonlinear strength. The time 
averaged probability of the exciton at site n is defined 
as 



< P n >= lim 



\C n \ 2 dt. 



(8) 



Therefore, we will look at the quantity |C„| 2 or < P n > 
or both of them for various situations to examine the 
occurrence of the self-trapping transition. 



III. RESULTS AND DISCUSSIONS 

First of all, we discuss the results for the Cayley tree 
with a single nonlinear impurity. The initial condition 
is set at the impurity site (zeroth site in Fig. 1). 




FIG. 3. The probability of the particle to be at the impu- 
rity site (Po (t)) as function of time for various values of non- 
linear parameter is shown. The long dashed curve, dashed 
curve, solid curve, dotted curve and the dotted-dashed 
curve correspond, respectively, to x = 4.75, 4.76, 4.77, 4.8 
and 5.0. 



The probability of finding a particle at the impu- 
rity site (initially populated site) is obtained by solv- 
ing Eq. (2), and the results are plotted in Fig. 3 as a 
function of time. Different curves correspond to the 
different values of the nonlinear parameter \- The 
long dashed, dashed, solid, dotted, and dotted-dashed 
curves correspond to the nonlinear strength of \ = 
4.75, 4.76, 4.77, 4.8 and 5.0, respectively. The con- 
nectivity for the Cayley tree considered here is K = 2. 

It is observed from Fig. 3 that, for \ = 4.75 and 4.76 
(the long-dashed and dashed curves), the probability 
of finding the particle at the impurity site decreases 
rapidly and then approaches to zero as time increases. 
This implies that the particle goes away from the im- 
purity site, i.e., the particle becomes fully delocalized. 
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However, the situation is drastically different for \ — 
4.77 (solid curve). The probability of finding the par- 
ticle at the initially populated site decreases down to 
about 0.35, then increases up to 0.61, and oscillates af- 
terward between 0.35 to 0.61. Thus, the average prob- 
ability of finding the particle at the initially populated 
site is approximately 0.48. For higher value of x, the 
probability of finding the particle at the initially pop- 
ulated site increases as is obvious from the dotted and 
dotted-dashed curves in Fig. 3. Therefore, we observe 
that there is a distinct critical value of x near (or be- 
low) 4.77, below which the particle escapes from the 
initially populated site and becomes fully delocalized, 
while above which the particle is most likely trapped 
at the initially populated site. 
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with a ~ 3.41 and (3 ~ 0.484. It should be noted that 
the critical value Xcr — 3.41 for K = 1 is prominently 
different from the critical value for a one-dimensional 
system. This implies that the power-law is valid only 
for K > 2. We therefore conclude that the geometry 
of the host lattice results in different critical values of 
X for different values of K > 2, obeying the power-law 
behavior in Eq. (9). 
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FIG. 5. The critical values of \ f° r self-trapping tran- 
sition in a Cayley tree with a single nonlinear impurity is 
plotted as a function of the connectivity of the Cayley tree. 
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FIG. 4. The time averaged probability of finding the par- 
ticle at the impurity site (i.e., < Po >) of the Cayley tree 
of K = 2 plotted against x- Sharp transition at \ = 4.77 
is found. 



The time averaged probability is also plotted in 
Fig. 4 as a function of the nonlinear parameter x- The 
sharp transition of the < Pq > (the time averaged prob- 
ability at the impurity site) at x — 4.77 is also clear 
from the Fig. 4. 

The self-trapping transition for the K = 1 case, i.e., 
for a single nonlinear impurity embedded in a one- 
dimensional lattice has been studied in detail by Dun- 
lap et al. pH , and the transition is found at x — 3.205. 
However, we notice that the self-trapping transition 
for a sinlge impurity embedded in the Cayley tree with 
if = 2 is sharper and clearer when compared with the 
case for the one-dimensional chain. 

In order to observe the dependence of the critical 
value of x a s a function of the connectivity K of the 
Cayley tree, we plotted in Fig. 5 Xcr for various val- 
ues of K > 2 on a double logarithmic scale. The data 
points lie surprisingly well on the straight line, imply- 
ing that the critical value obeys the power-law 

Xcr = a K\ (9) 
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FIG. 6. The time averaged probability of finding the pat- 
icle at the zeroth impurity site on a Cayley tree of K — 2 
with a dimeric nonlinear impurity plotted as a function of 
the nonlinearity parameter \- Sharp transition is observed 
at x = 4.95 



We now consider the dimeric impurity embedded in 
the Cayley tree with connectivity K = 2. The parti- 
cle is initially populated at one of the impurity sites. 
The probability of finding the particle at the initially 
populated site is calculated by solving the Eq. (7) for 
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various values of and the time averaged result is 
plotted in Fig. 6. The sharp transition is found at 
X cr = 4.95. One interesting observation found from 
Fig. 6 is that there is no precursor (peak) in the < Pq > 
before the permanent self-trapping transition occurs at 
X — 4.95, unlike the case in the one-dimensional sys- 
tem with two impurities (Fig. 5 in Ref. jp3[ ) for which 
a peak is observed at \ — 3.2 just before the perma- 
nent transition occurs at \ — 4.22. Thus, the par- 
ticle in the Cayley tree is always influenced by both 
the impurities present in the host whereas the parti- 
cle does not feel the presence of the second impurity 
in the one-dimensional system below x — 3-2. ^From 
Eq. (7) we note that the Cayley tree of connectivity 
K with dimeric impurity reduces to a one-dimensional 
chain with the hopping element between the impurity 
sites reduced from untiy to V c s = . For the Cayley 
tree K > 2, the hopping element between the impu- 
rity sites in the transformed one-dimensional system 
becomes less than or equal to Since the peak be- 
fore the permanent transition disappears in the case 
of Cayley tree, we suspect that there must be a crit- 
ical value for the hopping clement (say, V£g) between 
the dimeric impurity sites of a one-dimensional chain 
while keeping the other hopping elements unity such 
that the peak disappears for V c s < V^g. In order to 
verify this, we perform the numerical calculation of the 
time averaged probability for the particle to be at the 
initially populated impurity site of a one-dimensional 
chain with a dimeric nonlinear impurity for various val- 
ues of V e s (the hopping element between the impurity 
sites) while keeping the other elements unity. We inter- 
estingly observed that the gap (Ax) between the peak 
and the permanent transition point reduces as the V e s 
decreases and eventually vanishes at a critical value of 
V^ff = ^ (as shown in Fig. 7), i.e., the peak before 

the permanent transition dieappears when the system 
is equivalent to a Cayley tree with connectivity equal 
to 2. This implies that the structure of Cayley tree 
is responsible for the absence of the peak before the 
permanent transition. 
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FIG. 7. Plot of the gap Ax between the peak and the 
permanent transition point in the time averaged probability 
of the particle at the initially populated impurity site in a 
one-dimensional chain with a dimeric nonlinear impurity 
against the inverse of the square of the hopping element 
between the impurity sites. Other hopping elements are 
set to be unity. 



The critical value of x f°r various values of the con- 
nectivity of the Cayley tree is calculated and found to 
obey the same power-law behavior as in Eq. (9), but 
with different values of a and (3, i.e., a ~ 3.554 and 
(3 ~ 0.465. We note that the critical value of x f° r t ne 
self-trapping increases due to the presence of one more 
impurity. For a one-dimensional system the critical 
value is found to be 4.22 which deviates prominently 
from the value estimated from Eq. (9) with the values 
of a and (3 mentioned above and K = 1. This is again 
due to the different geometry of the host lattice. 

We now consider the perfectly nonlinear Cayley tree. 
The particle is initially placed at any arbitrary site. 
The sharp self-trapping transition is also observed in 
this case for various values of K. The critical value 
of the nonlinear parameter Xcr for self-trapping tran- 
sition again obeys the same power-law behavior as in 
Eq. (9) but with a ~ 3.766 and (3 ~ 0.445. Here we 
find that the critical values are larger than those for the 
case of a single and dimeric impurity; however, the dif- 
ference is not appreciably large. Therefore, it appears 
that the self-trapping transition is influenced by only 
few nonlinear neighbors around the initially populated 
impurity site. 

Finally it is worthy to mention that the formation 
of stationary localized state in the Cayley tree due to 
a single nonlinear impurity, dimeric nonlinear impu- 
rity and also in a perfectly nonlinear Cayley tree has 
been studied by Gupta and Kundu [p5| , p4| earlier and 
in the above three cases the bifurcations in the sta- 
tionary states were observed which in turn support the 
occurrence of the self-trapping transitions in the work. 

IV. SUMMARY 

The self-trapping transition due to a single and a 
dimeric nonlinear impurity embedded in the Cayley 
tree is studied. Furthermore, the self-trapping tran- 
sition in a perfectly nonlinear Cayley tree is observed. 
Very sharp self-trapping transition is observed for all 
systems considered here. The geometry of the host 
lattice is responsible for such a sharp transition. The 
critical value of the self-trapping transition increases as 
the number of nonlinear impurities in the host lattice 
increases. The critical value of x is found to obey a 
power law against the connectivity of the Cayley tree 
for all cases. Results are compared with those for the 
one-dimensional system. 
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